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Abstract 

In this work we present an atlas of composite spectra of 
galaxies based on the data of the Sloan Digital Sky Sur- 
vey Data Release 7 (SDSS DR7). Galaxies are classified 
by colour, nuclear activity and star-formation activity 
to calculate average spectra of high signal-to-noise ra- 
tio and resolution (S/N ~ 132 - 4760 at AA = 1 A), 
using an algorithm that is robust against outliers. Be- 
sides composite spectra, we also compute the first five 
principal components of the distributions in each galaxy 
class to characterize the nature of variations of individ- 
ual spectra around the averages. The continua of the 
composite spectra are fitted with BC03 stellar popula- 
tion synthesis models to extend the wavelength coverage 
beyond the coverage of the SDSS spectrographs. Com- 
mon derived parameters of the composites are also calcu- 
lated: integrated colours in the most popular filter sys- 
tems, line strength measurements, and continuum ab- 
sorption indices (including Lick indices). These derived 
parameters are compared with the distributions of pa- 
rameters of individual galaxies and it is shown on many 
examples that the composites of the atlas cover much 
of the parameter space spanned by SDSS galaxies. By 
co-adding thousands of spectra, a total integration time 
of several months can be reached, which results in ex- 
tremely low noise composites. The variations in red- 
shift not only allow for extending the spectral coverage 
bluewards to the original wavelength limit of the SDSS 
spectrographs, but also make higher spectral resolution 
achievable. The composite spectrum atlas is available 
online at |http : / / www . vo . elte . hu/ compositeatlas 
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1 Introduction 

With the advance of large scale spectroscopic surveys, 
spectrum stacking has become a widely-used technique 
in extragalactic astrophysics to derive composite spectra 
of high signal-to-noise ratios in order to enable the anal- 
ysis of spectral features that are washed away by noise 
in single object observations. Recent work have concen- 
trated on averaging and analysing single classes of ce- 
lestial objects, for example luminous red galaxies (22[, 
high-redshift galaxies galaxies in clusters [2(J, star- 
forming galaxies with varying inclination |48| , or quasars 

As the spectroscopic survey of the SDSS Legacy Sur- 
veys has finished, it is time to summarize the collected 
information. In this work, we focus on galaxies of vari- 
ous types, and present a comprehensive atlas of robustly- 
averaged composite spectra and spectral principal com- 
ponents of galaxies classified by many criteria. We use 
the original observed spectra as they were reduced and 
published in SDSS DR7 [J. 

1.1 Overview of the classification 

Galaxies are automatically classified following a revised 
classification scheme based on the work of Lee et al. [IH . 
In their work, galaxies were classified by colour, nuclear 
activity, star-formation activity and morphology. 

We separate galaxies into three colour classes: red, 
green and blue, green refe ring to green valley galaxies as 
defined by Strateva et al. [40|. Furthermore, galaxies are 
classified by star formation activity and nuclear activity 
into six classes: passive, measured Ha emission, star- 
forming, LINER and Seyfert galaxies; the sixth category 
is active galaxies with star-formation signatures hereafter 
denoted as AGN + Hn. Star-forming galaxies, LINERs 
and Seyferts are classified using the Baldwin-Phillips- 
Terlevich (BPT) diagnostic diagram [j| , and the segrega- 
tion curves determined by Ke wley et al. [3l[ , Kauffmann 
et al. [H| and Kewley et al. |32| . We omit morpholog- 
ical classification of the galaxies for reasons detailed in 
Sec. G2 where we also describe the classification scheme 
in detail. 
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1.2 Robust average and principal com- 
ponents 

To compute the average spectra, we use the robust sta- 
tistical method of Budavari ct al. [TH)]. This method not 
only gives the robust average and variance of a large, 
streaming set of data vectors (the spectra, in our case), 
but gives the first few principal components as well. Al- 
though we are primarily interested in the average spectra, 
the eigenspectra can help characterise the variations of 
objects within the subsets of the sample. 

Principal Component Analysis (PCA or Karhunen- 
Loeve Transformation) has been widely used in astro- 
nomical spectrum analysis in the last two decades. The 
basic idea of PCA is to derive an orthogonal basis from 
multidimensional data where the basis vectors point to 
the directions of maximum variance of the original data 
vectors. It was demonstrated that PCA can be success- 
fully applied to high-dimensional optical and UV galaxy 
spectra for classification purposes [17[. Real-life scenar- 
ios, when dealing with noisy and gappy data where cer- 
tain spectral bins are masked out due to bad observations 
or other reasons, require more advanced techniques 
An application of this robust technique to SDSS spec- 
tra was performed for galaxies and quasars by Yip et al. 
461 147| . A novel robust algorithm was developed by Bu- 
davari et al. flij which combines the robustness and gap- 
handling capability of the former methods with the abil- 
ity to handle extraordinarily large datasets by perform- 
ing the calculations iteratively. This new method signif- 
icantly reduces the computational and memory require- 
ments of the calculations and enables handling higher 
dimensional data (thus processing higher resolution spec- 
tra). 

1.3 Use cases of the spectrum atlas 

We would like to emphasize two use cases of this spec- 
trum atlas. One of the keys to the success of template- 
based photometric redshift estimators is the good set of 
spectral templates used to fit the broadband spectral en- 
ergy distributions of observed galaxies. Spectral tem- 
plates from models are often used, but they can miss 
some important features of real spectra. Also, models 
have to be highly consistent with the observations in 
terms of integrated colours. As an alternative to mod- 
els, individual empirical spectra are used, but it is not 
easy to get a high S/N set of spectra that represents all 
types of galaxies. Our composites can provide a high 
S/N, reliable and representative template set for photo-z 
applications. 



The other important use-case of our library is gener- 
ating mock catalogues. We expand the individual galaxy 
spectra on the eigenbases of of the various galaxy classes 
to determine the so called cigcncocfficients. The best fit- 
ting Gaussian curves to the distribution of these eigen- 
cocfficients are determined and the parameters of the 
Gaussian fits (m, a) are published, along with the eigen- 
spectra. By using these parameters of the distributions, 
linear combinations of the average composite spectra and 
eigenspectra can be easily constructed to further extend 
the parameter space coverage of the atlas. 

1.4 Structure of the paper 

In Sec. [2] we present the data we used for the analysis. 
In Sec. |3] we describe our classification scheme in detail. 
We present the details of the averaging method and the 
composite spectra in Sec. [4] where we also discuss the 
performance and robustness of the averaging algorithm. 
In Sec. [5] we fit stellar population synthesis models to 
the composites and present colours, emission line mea- 
surements and Lick indices of them. By comparing the 
derived physical parameters with the parameter distribu- 
tions of individual galaxies, we show that the composites 
of the atlas are representative to the whole SDSS galaxy 
ensemble. We summarize our work in Sec. El 



2 The data 

We analyse galaxy spectra from SDSS DR7 [l[ as reduced 
by the spectro pipeline [3!|, version 2.6. All spectra (in- 
cluding duplicate observations) are selected for analysis 
which were identified as galaxies by the spectro pipeline 
and have more than 3800 well-measured spectral bins. 

Photometric data are taken from the SDSS DR7 "best" 
dataset as measured by the SDSS photo image reduction 
pipeline (39j . 

K-corrections and absolute magnitudes are taken from 
the New York University Value-Added Galaxy Catalogue 
Q. K-corrections are calculated for z = 0.1, absolute 
magnitudes are based on the concordance ACDM cos- 
mological model with parameters of Hq = 100 km s _1 , 
Q. M = 0.3 and n A = 0.7. 

For line measurements of observed spectra, we use the 
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1 We apply the term "streaming" to data that is too large to fit 
into the memory so iterative loading from the background storage 
is necessary. 



MPA/JHU Value-Added Galaxy Catalogue 

Where not otherwise noted, galaxy luminosities are 
quoted in Petrosian magnitudes. Wavelength values are 
expressed in vacuum throughout the paper. 

3 Classification scheme 

We base our galaxy classification scheme on the work 
of Lee et al. [33j , but define several additional classes as 
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well. In the their scheme, galaxies were classified by mor- 
phology (early- type and late- type) , colour (red and blue) 
and nuclear activity (passive, star-forming, LINER and 
Seyfert). We extend this scheme by classifying galaxies 
by the strength of the Ha emission (completely passive 
galaxies and those showing some nebular emission but 
not enough to classify them using the BPT diagram) . We 
also introduce a transitional class between star-forming 
galaxies and active galactic nuclei (AGNs) based on the 
works by Kauffmann et al. [28| , Kewley et al. [3ll [32| . 
Kauffmann et al. [28[ termed this latter class 'compos- 
ite'; we will refer to these galaxies as 'AGN + Hn' to 
avoid confusion with composite spectra. The colour clas- 
sification is extended to treat green-valley galaxies sepa- 
rately. Volume-limited cuts are imposed on all subsam- 
ples as described in Sec. 13.31 For the most numerous 
galaxy classes, we also define further classes based on 
colour cuts in Sec. 13.41 and term these classes 'refined 
star-forming' or 'SF', and 'refined red' or 'RED'. The 
number of galaxies used to calculate the composites, the 
wavelength coverage and the estimated signal-to noise 
ratio are summarized in Tab. [T] 

In contrast to Lee et al. [331 ] . we omit morphological 
classification of the galaxies for the following reasons. 
Morphological classification relies on the radial colour 
gradient of galaxies which was not measured perfectly 
by the SDSS photo pipeline for the whole spectro sam- 
ple. Thus, galaxy classes defined by morphology would 
contain about 60 per cent fewer galaxies: only those 
with good radial colour gradient measurements. Using 
morphological information would reduce the number of 
the classified galaxies to a number that would not al- 
low reaching the desired low noise level and resolution. 
Also, the size and magnitude dependent incompleteness 
of the morphological classification would introduce un- 
wanted bias into the data set. 

We further explain the classification criteria in detail 
in the following subsections. 

3.1 Classification by colour 

We organise the galaxies into three colour classes based 
on their redshifts and g ~ r colours. The three colour 
classes are the following: 

• Red galaxies: Having primarily old stellar popula- 
tions or galaxies with significant intrinsic extinction. 

• Green galaxies: Those belonging to the green val- 
ley, thought to form the transitional region between 
active star- formation and passivity 0, 0, |4(j ■ 

• Blue galaxies: Mostly actively star-forming galaxies. 

Fig. Q] displays the colour classes we define. Galaxies 
above the solid red curve are classified as 'red', galaxies 
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red A+Hn 
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483 


red S 


3350 
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SF 3 
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SF 4 


3550 
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4436 
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SF 5 


3550 
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Table 1: Basic parameters of the composites by activity 
and colour class. 



3 




I 



0.6 



0.4 I 



0.2 



0.00 0.05 0.10 



0.15 0.20 
rcdshift 



0.25 0.30 0.35 



Figure 1: Definitions of the colour cuts. Density plot: 
logarithm of the distribution of the galaxy sample, nor- 
malised in each redshift bin. Red galaxies: above the 
solid red line; blue galaxies: below the dashed blue line; 
green valley galaxies: between the two coloured lines. 
The data points and the dashed white curve represent 
the peak of the red sequence and our best fit to them 
as a function of redshift. Sec online edition for colour 
version of this plot. 



below the blue curve are classified as 'blue' and green 
galaxies lie inbetween the two. These curves were deter- 
mined by fitting Gaussians to the distribution of galax- 
ies by the g — r colour in bins of Az = 0.05 between 
0.0 < z < 0.35. The dashed white line is the center of 
the red sequence, the solid red line is the 1.2<r envelope of 
it, and the blue line is at 3.5<r from the red sequence. To 
construct this plot, all galaxies were used independently 
of their apparent or absolute magnitudes. 

In Fig. [1] colour curves are defined by the a param- 
eter of the fitted Gaussians as explained above. Since 
the photometric error increases with increasing redshift, 
it makes the scatter in the red sequence higher. This 
causes the colour curves to turn over around z ~ 0.2 
and be more restrictive for green-valley and blue galax- 
ies, allowing more variance in the colour of red galaxies. 
Because of the relatively low number of observed green- 
valley and blue galaxies at these redshifts, this would 
have no effect on the average composites, but more high- 
redshift red galaxies could be used for the analysis. In 
the present analysis we set the redshift limits of the sub- 
sample volumes below z < 0.17. 

3.2 Classification by activity 

Our classification scheme by star-formation and nuclear 
activity is based on the works by Kewley et al. (3l| and 



Kauffmann et al. 
Kewley et al. [32 
based on activity: 



|. LINER/Seyfert separation follows 
We define the following six classes 



• Passive galaxies: No measurable Ha emission. 

• Ha measured: Hydrogen emission is measured, but 
other lines are not strong enough for classification. 

• Star-forming galaxies: with clear star-formation sig- 
natures but no active nucleus. 

• AGN + Hn galaxies: with AGN and extreme star- 
burst signatures at the same time. 

• LINER galaxies: pure LINER signatures pHj . 

• Seyfert galaxies: pure Seyfert signatures. 



In Fig. [5] we plot the BPT diagram of the galaxies in 
our sample. The segregation line (solid red curve) be- 
tween star-forming galaxies (below the curve) and AGN 
(above the curve) is: 
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The limit of extreme starburst (blue dashed line) is 
determined as: 
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Figure 2: BPT diagram for our sample of SDSS galax- 
ies. The density plot is the square root of the number 
density of SDSS galaxies with high enough line S/N mea- 
surements (see text). The solid red curve indicates the 
segregation line between star forming galaxies and active 
galactic nuclei. The blue dashed curve is the extreme 
starburst line. The dashed rectangle marks the same 
area as displayed on Fig. [T5] See text for references. 
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Figure 3: Seyfert/LINER separation. The density plot 
is the square root of the number density of our sample of 
SDSS galaxies classified as AGN (without star-formation 
signature) using Fig. [2] The line separates the two sides 
of the bimodal distribution. See text for references. 



Although their number is negligible, we excluded 
galaxies with extreme values of line ratios, i.e those that 
are outside the limits of the plot in Fig. [2] 
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The latter line separates pure AGN (above the curve) 
from AGN showing star-formation signatures (between 
the two curves). 

Fig. |3] shows the distribution of galaxies classified as 
AGN (without star- formation signatures). We apply the 
following cut to distinguish LINERs (below the line) from 
Seyfert galaxies (above the line): 



log 



10 



[ Qui] 

H/3 



= 1.18 log 10 ^1 + 1.4 



Ha 



Only galaxies with Ha measured at S/N > 3 are clas- 
sified as star-forming, while only those with all Ha, H/3, 
[Ol], [Oni] and [Nil] measured at S/N > 3 arc classified 
as AGN (either AGN + Hn, LINER or Seyfert). For 
those galaxies that do not have strong enough emission 
lines to make it onto the BPT, we define two further 
classes: completely passive with no detectable emission 
lines and passive with measurable Ha. 



3.3 Volume-limited samples 

In order to avoid the Malmquist bias, we define redshift 
and absolute magnitude ranges for every galaxy class, 
trying to maximize the number counts of the volume lim- 
ited samples. 

Although setting restrictive constraints on the redshift 
and absolute magnitude reduces the wavelength coverage 
of the composites, it is necessary to get physically mean- 
ingful results. The spectroscopic targeting algorithms of 
SDSS 21, H, 41 1 set the sharp faint-end r-band magni- 
tude limit at 17.7 mag for the main galaxy sample. On 
the other hand, bright-end incompleteness is much less 
well controlled. One of the main reasons for this incom- 
pleteness is that nearby galaxies are too big in apparent 
diameter compared to the spectroscopic fibre diameter. 
We estimate the bright-end apparent r-band magnitude 
limit to be around of 13.5 when determining the volumes. 

We set the volumes as plotted in Fig. |4] and summa- 
rized in Tab. d 
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Figure 4: Definitions of redshift and r-band absolute magnitude M r cuts for every galaxy activity class. Density plot: 
square root of the galaxy density within the activity class; solid red line: faint-end magnitude limit corresponding 
to 17.77 mag, gray line: estimated bright-end magnitude limit corresponding to about 13.5 mag. The cuts were 
manually set to contain as many galaxies as possible considering the magnitude limits of the SDSS spectro sample. 
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passive 


0.03 


0.14 


-20.5 


-21.5 


Ha measured 


0.03 


0.14 


-20.5 


-21.5 


star-forming 


0.03 


0.08 


-19.2 


-21.0 


AGN + Hn 


0.03 


0.14 


-20.5 


-21.5 


LINER 


0.04 


0.17 


-21.0 


-22.0 


Seyfert 


0.03 


0.14 


-20.5 


-21.5 



a 



Table 2: Rcdshift and absolute magnitude ranges for 
each galaxy activity class. 



3.3.1 Effect of incomplete sampling 

We note that the continua of average composite spec- 
tra are physical in the sense that they represent a non- 
negative linear combination of real stellar populations. 
This makes it meaningful to fit composites with theoret- 
ical population synthesis models. However, incomplete 
sampling of galaxies for composite calculation has the 
effect of making the composites non-physical. 

While we consider spectral evolution negligible in the 
small redshift intervals we choose, care must be taken 
about the dependence of spectra on the absolute mag- 
nitude of the objects. For instance, choosing a lower 
redshift cut might enable extending the wavelength cov- 
erage of the composites redwards, but bright-end incom- 
pleteness would lead to the dominance of lower luminos- 
ity galaxies which (in case of blue galaxies, for example) 
tend to be bluer. This would result in a downturn of the 
composite spectrum as it can be seen in Fig. [SJ 

3.4 Refined classification of blue and red 
galaxies 

Template-based photometric redshift (photo-z) estima- 
tion algorithms (for a recent review, see Hildcbrandt 



et al. [26[ and references therein) rely on good qual- 
ity empirical galaxy spectra or theoretical models to fit 
the broadband observations with; though certain tech- 
niques exist to improve the reliability of the models 14 1 . 
Most often, the problem with using theoretical models 
for photo-z is the mismatch between the colours calcu- 
lated from the models and the colours of the observed 
galaxies. Also, photo-z techniques are usually based on 
a small number of templates in order to reduce compu- 
tational complexity and to avoid spurious results due to 
fitting too many free parameters of the not-so-reliablc 
models. 

The most numerous galaxies (passive red ellipticals 
and star-forming blue spirals) show significant variations 
in their colours. Using only a few templates might cause 
the photo-z estimates to be coarse because the templates 
are unable to characterize this intrinsic colour variance of 




Figure 5: The effect of incomplete sampling. Using low 
rcdshift spectra would allow extending the composites 
redwards but bright-end incompleteness causes the aver- 
aged spectrum to turn down at longer wavelength making 
the composite non-physical. Solid black line: composite 
spectrum calculated from an incomplete sample, dashed 
red line: best fit model. (Continuum fitting was done in 
the 3500 - 9000 A range). 



the objects. More reliable empirical templates, however, 
could be used to refine photometric rcdshift estimation. 
In order to address this issue, we further refine the most 
numerous galaxy classes by colour and calculate their 
composite spectra. 

3.4.1 Blue galaxies 

In the case of star- forming galaxies (mostly spirals and ir- 
regulars) , the intrinsic colour variance is caused by many 
effects. Star- formation rate is primarily responsible for 
the restframe colours of galaxies, while dus t content and 
inclination play an important role too 



481 ] . Also, since 



fibre spectrographs with limited aperture only observe 
the central, redder region of the galaxies, the distance of 
the galaxies and their bulge-to-disc ratio is an important 
factor causing colour variations in small aperture magni- 
tudes. We will denote these classes as 'SF' followed by a 
number. 

3.4.2 Red galaxies 

In the case of red galaxies (mostly ellipticals) , the origin 
of the colour variation is mainly an evolutionary effect 
caused by their passively evolving old stellar populations 
(for a comprehensive review, see Renzini [37j and refer- 
ences therein); the effect of dust content on their optical 
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colours is thought to be minimal. 

For the purposes of refined colour classification, we will 
use all galaxies of the former 'passive' and 'Ha measured' 
classes and denote these class as 'RED' followed by a 
number. 



3.4.3 The refined colour classes 

We define the refined colour classes based on the g — r 
colour distributions of star-forming and passive galaxies 
as plotted in Fig. El The histograms of the distribu- 
tions are fitted with Gaussians and the colour limits are 
defined at 0.5a and 1.5a distances from the mean. Alto- 
gether, five additional colour classes are defined, as listed 
in the following table: 



# 


limits 


1 




g — r < m - 


- 1.5a 




2 


777, — 


1.5a < g — r 


< 771 — 


0.5a 


3 


m — 


0.5a < g — r 


< 771 -+■ 


0.5a 


4 


m + 


0.5a < g — r 


< m 4 


1.5a 


5 




m + 1.5a < 


g — r 





Parameters m and a denote the mean and standard 
deviation of the fitted distributions, respectively. The 
fitted parameters are summarised in the following table: 





777, 


a 


SF 
RED 


0.48 
0.95 


0.16 
0.06 



3.5 Redshift distribution 

The redshift distribution of the galaxies is plotted in 
Fig. [7] for every galaxy activity class in redshift bins of 
Az ~ 0.01. It is clear from the curves that all sam- 
ples except red galaxies are flux-limited according to the 
spectroscopic targeting algorithm of SDSS in the investi- 
gated redshift range. The bright-end red galaxy sample 
is volume-limited to z = 0.35 [2l|- Fig. [8] shows the 
same redshift distributions after applying the absolute 
magnitude cuts to get the volume-limited samples. Note 
that the requirement to select volume limited samples 
does not allow maximizing galaxy counts by freely choos- 
ing the used redshift intervals; absolute magnitude limits 
also have to be taken into account. 




-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 



9.0 




g-r [mag] 



Figure 6: g—r colour distribution of star-forming galaxies 
(top panel) and passive galaxies (bottom panel). The 
red histograms show the measured distributions; black 
curves indicate the fitted Gaussians. The vertical dashed 
lines mark the colour cuts we defined at 0.5cr and 1.5<7 
distances from the mean. 



4 Composite spectra 
4.1 Preprocessing 

We selected spectra based on the criteria defined in Sec. [3] 
and preprocessed them before calculating the compos- 
ites. First, all spectra were corrected for foreground 
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Figure 7: Rcdshift distribution of the galaxies in each activity class before the absolute magnitude cuts. Colours 
represent the colour classes: red - red galaxies, green - green valley galaxies, blue - blue galaxies, and black - all 
galaxies regardless of their colours. The dashed vertical lines represent the redshift cuts, which are different for 
every activity class. 
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Figure 8: Redshift distribution of the galaxies in each activity class after the absolute magnitude cuts. Colours 
represent the colour classes: red - red galaxies, green - green valley galaxies, blue - blue galaxies, and black - all 
galaxies regardless of their colours. The dashed vertical lines represent the redshift cuts, which are different for 
each activity class. 
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galactic extinction using the dust map of Schlegel et al. 



35j. Sec- 



38[ and the extinction curve from O'Donnell 
ond, spectra were normalised as described in Sec. 14.21 
Third, spectra were converted to rest-frame and re- 
binned to a grid with the intended resolution (AA = 1.0, 
0.5 or 0.25 A). Q 

4.2 Normalisation 



To derive meaningful composites, one has to normalise 
the galaxy spectra in a robust way. Normalising spectra 
at a given wavelength is not appropriate for noisy mea- 
surements, thus we use the following technique. We cal- 
culate the median flux in the 4200 - 4300 A, 4600 - 4800 
A, 5400-5500 A, and 5600-5800 A wavelength intervals 
and use this value to scale spectra together. The wave- 
length intervals are chosen to be devoid of any strong 
emission lines. Spectra are re-scaled to have the median 
value of the flux in these intervals equal to one. 

We do not weight or otherwise scale the flux of the 
spectra based on the apparent or intrinsic luminosity 
of the objects or S/N of the observations. Although 
this could affect the signal to noise ratio of the com- 
posites negatively, the robust algorithm can successfully 
overcome this problem by automatically assigning lower 
weights to spectra with lower signal to noise values. 



4.3 Robust algorithm 

In order to take the significant noise and the varying sig- 
nal to noise ratio of the original data set into account, 
we use the robust averaging method of Budavari et al. 
[l5| to determine the average and variance spectra of the 
different galaxy classes. This algorithm not only com- 
putes the average and the variance iteratively by taking 
care of gaps in the measurements, but also yields the 
first few principal components of the spectrum distribu- 
tion (see Sec. I4.5[) . We will detail the efficiency of the 
robust method with respect to ordinary average and me- 
dian composite spectra in Sec. 14.61 

The composites calculated by the robust algorithm are 
plotted in Figs. IMTSIfor each activity and colour class. 

4.4 High resolution 

Because composites are combined from many spectra at 
varying redshifts, it is possible to get a couple of factors 
higher resolution than the resolution of the original mea- 
surements. The average error in redshift measurement 
of the high signal-to-noise (S/N > 30) galaxy spectra of 
SDSS is z elI = 1.6 x 10~ 4 , which allows room for over- 
sampling the composites. On the other hand, oversam- 
pling can introduce noise compared to the almost noise- 
less composites at 1 A binning. This noise can become 
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Figure 13: The K line of a galaxy composite spec- 
trum calculated at resolutions AA = 1,0.5,0.25 and 0.1 
A(from bottom to top). Higher resolution curves are 
smoother compared to the AA = 1 A curve. 



significant for sparsely populated galaxy classes such as 
the green valley galaxies. 

Another important factor that counteracts the efforts 
to get higher resolution is the intrinsic broadening of 
the spectral features of the individual spectra due to the 
high velocity kinetics in the galaxies. To overcome this 
problem, deconvolution of the Dopplcr broadening kernel 
from the original measurements would be necessary. We 
do not address this issue in the current work. 

Fig. [TBI shows the comparison of some spectral features 
of a galaxy composite with different oversampling (at the 
resolution of 1 A, 0.5 A, 0.25 A, and 0.1 A). Note how 
the line features become smoother with higher resolution. 
It is important, though, that oversampled spectra might 
not contain more information because the wide velocity 
dispersion washes out fine details. 

Evidently, linear interpolation could not be used to get 
the higher resolution composites form lower resolution 
ones. Higher-order curves could be used for oversampling 
low resolution composites to get smooth line profiles, but 
they are likely to introduce artefacts. Our method of 
oversampling overcomes this problem. 

4.5 Generating empirical mock cata- 
logues using eigenspectra 

Besides the robust average spectrum of every galaxy 
class, we also calculate the variance spectra. These can 
be used to determine the wavelength ranges with out- 
standing variance (usually the emission lines) and quan- 
tify the noise in each wavelength bin. The variance spec- 
tra are available on-line. However, if one is interested 
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Figure 9: Composite spectra of the different activity and colour classes. Activity class is indicated in the top left corner of each panel, the spectra of 
the different colour classes are offset by 0.5 along the y axis for clarity, and indicated by colours from top to bottom (black: all colour, red, green and 
blue). See the online edition for a colour version of this plot. 
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Figure 10: See caption of Fig. [9] for detailed description. 
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Figure 11: Sec caption of Fig. [9] for detailed description. 
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Figure 12: Composite spectra of the refined colour classes as described in Sec. 13.41 The curves are colour-coded from blue (top) to red (bottom) based 
on the g — r colour of the galaxies. See the online edition for a colour version of this plot. 



in the overall variance of the spectra (in contrast to the 
bin-by-bin variance) of a galaxy class, it is not enough to 
use the variance spectrum where bins are treated inde- 
pendently. In order to quantify the variance among the 
spectra of a certain galaxy class over the entire wave- 
length range, and take the covariance of the bins into ac- 
count, the first few principal components (cigenspectra) 
have to be calculated. Fortunately, our robust averaging 
algorithm yields the eigenspectra as a "by-product" of 
the averaging. 

We calculate the first five eigenspectra of every galaxy 
class, and expand each measured spectrum on this basis. 
In most cases, the distribution of the expansion coeffi- 
cients is very close to a Gaussian distribution centred on 
zero. We determine the standard deviation cr i g of the 
distribution of the eigencoefficients in each galaxy class. 
Then one can calculate the following spectra. 



F ± _ pavg 



_|_ U a . pcig 



(1) 



where F^ vs and F^ 1S are the average and one of the 
eigenspectra, respectively. The coefficient k is a constant 
factor that can be varied to get slightly different spectra. 

As an illustration, we plot for the refined classes of 
red and blue galaxies in Fig. [14] In case of the red galax- 
ies, the first eigenspectrum clearly determines the colour 
of the galaxy; almost the entire variation in the overall 
slope of spectrum is picked up by the first principal com- 
ponent. In case of the star-forming galaxies, most of the 
variance is in the emission lines, so one has to use the 
second eigenvector to generate spectra with colour vari- 
ations. As Fig. [Ml shows, this technique can be success- 
fully used to generate composites of any colour between 
the averages of the predefined classes. One drawback of 
the method is that adding higher order eigenspectra to 
the averages introduces some noise. 

The importance of this technique is that it allows gen- 
erating mock spectral libraries. By drawing the value 
of k of Eq. [1] from a normal distribution, the F^ 1 spec- 
tra will be representative to the whole ensemble of SDSS 
galaxies. Eigenspectra of the various galaxy classes are 
published online along with the parameters (mean and 
a) of the distributions of the expansion coefficients of in- 
dividual observations on the eigenbases. By using these 
parameters, linear combinations of the average compos- 
ite spectra and eigenspectra can be easily calculated to 
further extend the parameter space coverage of the pre- 
sented atlas. 

4.6 Performance of the robust method 

To quantify the performance of the robust method, we 
compare the robustly averaged composites with normally 
averaged and median composites. Averaging is thought 
to result in composites that are physical in the sense 



that they are constructed as non-negative linear combi- 
nations of real stellar continua. Median, while it is robust 
against noise, differs from averaging because it treats the 
spectral bins independently and might yield non-physical 
continua. One disadvantage of averaging spectra it that 
emission line ratios are not guaranteed to be preserved, 
whereas Vanden Berk et al. |43[ reported that taking the 
median spectra of quasars did not alter the emission line 
ratios. We will revisit the question of emission lines of 
average composites in Sec. 15.31 The main advantage of 
robust averaging over median calculation is that the for- 
mer can be done iteratively, and there is no need to store 
and sort all the data in memory. 

We plot the difference between the robust average and 
the non-robust composites in Fig. [15] and [16] (middle pan- 
els). The robust composite and the number counts used 
to calculate the composites in each wavelength bin are 
also plotted for reference in the top and bottom panels, 
respectively. 

Interestingly, for red galaxies, the median runs much 
closer to the robust average than the normal average, 
while in case of blue galaxies the differences are not as 
significant as in case of the red galaxies. Significant dif- 
ferences between the robust and non-robust composites 
are visible at very short and very long wavelengths, and 
also at emission lines. The differences between the con- 
tinua are possibly due to the low number counts in those 
wavelength ranges. For low number counts, non-robust 
methods tend to become severely biased. The robust al- 
gorithm assigns lower weights to galaxies with extremely 
strong emission lines, which is considered to be the rea- 
son behind the large difference between the robust and 
non-robust methods in the emission line regions. 



5 Composites on the map 

In this section, we use standard techniques to estimate 
the stellar continua, emission line strengths, absorption 
line indices and physical parameters of the composites. 
We follow the methods described by Brinchmann et al 
Kauffmann et al 



42|. 



[11| . Kauffmann et al. [28|, |29| and Tremonti et al 
These techniques were developed and used by the the 
cited authors to compose the MPA/JHU value-added 
galaxy catalogue!! of SDSS DR7. 

In this paper, we do not intend to make any conclu- 
sions about physical parameters of the galaxy samples 
based on the parameters derived from the composites. 
Our goal is simply to show that the composites calculated 
using our method cover much of the parameter space of 
the individual galaxies. To achieve this, we plot several 
standard diagrams widely used in the literature. These 



diagrams show the distributions of the parameters of in- 



^http : //www.mpa-garching .mpg. de/SDSS/ 
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Figure 14: Variance in red (left panel) and blue star-forming (right panel) galaxies. The solid curves show the 
composites of galaxies classified according to the refined colour classes (redder to bluer) . The shaded areas represent 
the 0.5er (red galaxies) and la (blue galaxies) variations about the average. See the online edition for a colour version 
of this plot. 
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Figure 15: Comparision of the various averaging methods for red galaxies with measured Ha. We plot the spectrum 
in the top panel for reference. The middle panel shows the difference between the robust average and the normal 
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Figure 16: Same as Fig. [15] but for star-forming blue galaxies. 



dividual galaxies of the SDSS or MPA/JHU value-added 
catalogues overplotted with the parameters of the com- 
posites. 

5.1 Continuum model fits 

The wavelength coverage of our empirically-derived high 
signal to noise composite spectra spans from about 3350 
A to 9050 A. The actual coverage depends on the red- 
shift cuts of the volume-limited samples of the galaxy 
classes; higher rcdshift samples allow extending the cov- 
erage blueward. In order to extend this limited cover- 
age, we fit the composites with theoretical model spec- 
tra generated using the single stellar population (SSP) 
synthesis models of Bruzual & Chariot [13( . Three sets 
of models were generated with different metallicity val- 
ues of Z = 0.008, 0.02, and 0.05. Each set contains 10 
single burst models with ages of t = 0.005, 0.025, 0.100, 
0.290, 0.640, 0.900, 1.4, 2.5, 5.0 and 11.0 Gyr based on 
the initial mass function (IMF) of Chabrier 11611 and the 
evolutionary tracks 'Padoval994' 0, US 13- We use 
the non-negative least squares (NNLS) method to deter- 
mine the best linear combination of models of different 
age (but same metallicity) to fit the composite spectra. 
To estimate the effect of intrinsic extinction by dust, we 
use the simple dust model of Bruzual & Chariot [1J 
which is described by the following formula: 



rx{t) 



ry(A/5500 A)" 07 for t < 10 7 yr 
^7y (A/5500 A)-°- 7 for t > 10 7 yr 



We let the value of the optical depth ry vary while 
keeping the value of fi = 0.3 fixed. In this particular dust 
model, /i controls the ratio of the effect of extinction on 
stars older/younger than 10 Myr. We also keep the value 
of the velocity dispersion of the stars a v a free parameter. 
For each composite, we repeat the fitting for each model 
set, and the metallicity of the best fitting model set is 
used. 

To compute the synthetic magnitude of the compos- 
ites in broadband filters where the wavelength coverage 
of the composites is not wide enough, we extend the em- 
pirical spectra with the best fitting models. We also use 
the same continuum fits to measure the emission line 
equivalent widths of the composites. 

5.2 Colour-colour trajectories 

As we intend to use the composite atlas as a basis for 
template-based photometric redshift estimation, it is im- 
portant that the composites cover much of the colour 
space spanned by the SDSS galaxies. Fig. [T7] and [15] dis- 
play the SDSS g — r; r — i colour-colour trajectories of 
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Figure 18: Colour-colour diagrams of the composites for 
the SDSS g — r and r — i colour indices between redshifts 
< z < 0.7. The coloured curves represent the compos- 
ite spectrum of the different refined colour classes (top 
panel: red galaxies, bottom panel: star-forming galax- 
ies), as defined in Sec. 13.41 Dots along the curves mark 
redshifts z = 0.0, 0.1, etc. The greyscale map in the back- 
ground shows the distribution of the photometric colours 
of SDSS spectro galaxies. 



composites for every galaxy class we defined. The trajec- 
tories cover the redshift range of 0.0 < z < 0.7, and they 
nicely overlap with the measured galaxy distributions. 
This is a good confirmation that the robust averaging al- 
gorithm does not introduce any bias into galaxy colours. 

It is observable in Fig. [171 that the bluest galaxies (g — 
r < 0.5) are not covered by the composites of the main 
galaxy activity classes. The refined classification of star- 
forming galaxies however yields composites of extremely 
blue colours, as it can be seen in[T8l 



5.3 Emission line measurements 

As averaging spectra may have the effect of altering line 
ratios, we have to verify whether the emission lines of 
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Figure 19: BPT diagram of the composite spectra. The 
density plot is the logarithm of the number density of all 
SDSS galaxies. The solid red curve indicates the segrega- 
tion line between star forming galaxies and active galac- 
tic nuclei. The blue dashed curve is the extreme starburst 
line. The loci of the composites are indicated with points: 
rectangles - star-forming galaxies; diamonds - AGN+Hn 
galaxies, triangles - LINERs, circles - Sey forts. The dif- 
ferent colours denote the different colour classes; black 
points mark the loci of composites of all galaxies (re- 
gardless of their colours) in a given activity class. Colour 
crosses mark the loci of the star-forming galaxies, accord- 
ing to the refined classification. See the online edition for 
a colour version of this plot. 



the composites are physical or not. Fortunately, in our 
case the robust averaging algorithm does not alter the 
line ratios severely, as we will show below. 

First, emission line equivalent widths are determined 
by using the best fit continuum model as described in 
Sec. 15.11 Then, the BPT diagram is plotted for the com- 
posites of the star forming, AGN + Hn, LINER and 
Seyfert galaxies in Fig. Q1J] The loci of the composites 
overlap with the line ratios calculated from the individual 
galaxy measurements. It is important to mention, how- 
ever, that the variance in line ratios of galaxies is much 
larger than what a few composite spectra can character- 
ize. 

One interesting feature of Fig. [19] is that Seyferts are 
located along a well-defined line parametrized by the g—r 
colour: redder Seyferts tend to have higher [Oiii]/H/3 
ratios. This is consistent with the findings of Kewley et 
al. (32[: the stellar populations of AGNs are the older 
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Figure 17: Colour-colour diagram of the composites for the SDSS g — r and r — i colour indices between redshifts 
< z < 0.7. The coloured curves represent the composite spectra of the different colour classes we defined in 
Sec. 13.11 The black curve represents the composite spectra calculated from galaxies with any colour. Dots along 
the curves mark redshifts z = 0.0,0.1, etc. The greyscale map in the background shows the distribution of the 
extinction corrected photometric colours calculated from the fibre magnitudes of SDSS photo galaxies, which is 
significantly deeper than the spectro set. 
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the further away they are located from the star-forming 
line on the BPT diagram. For AGN + Hn galaxies, no 
similar behaviour is readily visible. 

In case of the star-forming galaxies, we also plot 
the results from the refined colour classification. A 
clear trend is observable here too: bluer galaxies have 
higher [Om]/H/3 ratios and lower [Nil] /Ha ratios. Bluer 
star- forming composites (crosses), however, tend to be 
slightly above the star-forming line of the MPA/JHU re- 
sults (density plot). The most trivial explanation to this 
effect is that averaging alters the line ratios and this bias 
is more significant in case of stronger emission lines. 

Emission line equivalent widths of all composites are 
listed in Tab. [3] -[5] and are available on-line. 



5.4 Continuum indices 



We calculate numerous absorption line indices of the 
composites: Lick indices 
Hanes, sec Huchra et al. 



44 |45| , BH indices (Brodie & 
27| and DTT indices 0. In- 



dices characterizing the size of the 4000 A break are from 



Bruzual G. [13 and Balogh et al. |6j. Absorption indices 



for all composites are summarized in Tab. and are 
available online. 

One widely-used diagram to separate star-forming, 
starburst and post-starburst galaxies from passive ones 
based on absorption line features is the plot of the 
D4000„ -US a plane [11,13^. As an illustration to the 
parameter space coverage of our composite atlas, we plot 
the D4000„ - R5 A diagram in Fig. [20] and EU Ac- 
cording to Kauffmann et al. 30 . in this plane, galaxies 
are located along a straight line: low-mass star-forming 
galaxies are close to the top left corner of the plot while 
high-mass passive galaxies are in the bottom right cor- 
ner. Intermediate-mass galaxies (10 10 < M Q < 10 11 " 5 ) 
show some bimodality and populate the two big clumps 
as well as the line between them. Starburst galaxies are 
located slightly above the straight line [29[ . The compos- 
ites clearly follow the relationship among colour, Ha and 
^ISa '■ passive red galaxies with no or weak Ha emissions 
are in the lower right corner, blue star-forming galaxies 
with strong Ha emissions are in the top left corner of the 
plots. Most AGNs fall into the intermediate-mass region, 
only the reddest LINERs overlapping with the clump of 
the completely passive red galaxies. In terms of recent 
star-formation activity, the two interesting galaxy classes 
are passive blue galaxies and blue LINERs. These galaxy 
classes are likely to contain many post-starburst galaxies 
as they are slightly above the linear distribution of the 
majority of galaxies. 



5.5 Physical parameters from Bayesian 
analysis 

The model fitting method described in Sec. 15 . II is known 
not to yield correct values for the physical parameters 
of the stellar populations building up the galaxies, and 
more sophisticated techniques are required to estimate 
the age and the metallicity of the galaxies. 

To estimate the physical parameters of the continua of 
the composites, we entirely adopt the method described 
in Sec. 4 of Kauffmann et al. 29]. They used a Bayesian 
approach to sample the parameter space of the BC03 stel- 
lar population synthesis models [l3| paying special atten- 
tion to generate physically plausible star formation his- 
tories and constrain other parameters (metallicity, dust 
content). Based on the stochastic library of star forma- 
tion histories, they generated a huge library of model 
spectra. The star formation histories were characterized 
by a continuous rate starting at iform and decreasing ex- 
ponentially with the time-scale 7 superimposed by ran- 
dom, short-term bursts. 

To estimate physical parameters of a single spectrum, 
models of the library are individually fitted to the com- 
posites and from the goodness of these fits a probability 
of the model parameter set is calculated. This method 
gives not just a single best-fitting model to the spectrum 
as the NNLS method described in Sec. 15.11 does, but 
yields the whole probability distribution of the model 
parameters and physical properties. As certain physi- 
cal parameters are very hard to constrain based on opti- 
cal spectroscopy, it is much better to characterize galaxy 
spectra with the probability distribution of their param- 
eters instead of a single best-fit model. The distribu- 
tions of the physical parameters of our composites are 
determined using the described method and the modes 
and medians of the distributions are listed in Tab. [9[ 
It is clear from the results that our composites cover a 
very broad range of stellar population ages, specific star 
formation rates and mass-to-light ratios. On the other 
hand, metallicities are weakly constrained and close to 
solar. For the same parameter measurements of individ- 
ual SDSS galaxies, see Brinchmann et al. 11 1. 



6 Summary 

We have presented a comprehensive atlas of composite 
spectra of SDSS galaxies. The signal-to-noise ratio of 
the composites can reach as high as S/N ~ 132 — 4760 
at AA = 1 A sampling. We have shown that the vol- 
ume limited sampling is very important to get physi- 
cally meaningful composites at very short and very long 
wavelength. We have argued that our robust averaging 
method is superior to both traditional averaging and me- 
dian calculation in terms of computational efficiency and 
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Table 3: Values of the emission line equivalent widths (in Angstroms) for every galaxy class, continued in Tab. 01 
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Table 4: Values of the emission line equivalent widths (in Angstroms) for every galaxy class, continued from Tab. |3l 
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Table 5: Values of the emission line equivalent widths (in Angstroms) for the refined galaxy classes. 
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Figure 20: Recent star-formation signatures in the composites. The density plot is the square root of the number 
density of all SDSS galaxies. The loci of the composites arc indicated with circles. The different colours denote 
the different colour classes; black points mark the loci of composites of all galaxies (regardless of their colours) in 
a given activity class. See the online edition for a colour version of this plot. 
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Figure 21: Recent star- formation signatures in the refined composites. The density plot is the square root of the 
number density of all SDSS galaxies. The crosses mark the loci of the composites based on the refined colour 
classification; see Sec. 13.41 
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Table 6: Values of the continuum absorption indices for every galaxy class, continued in Tab. [71 
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Na D 


x 

A 


2 


566 


2 


209 


1.282 


2.363 


2 


447 


2 


121 


2 


208 


3 


366 


3 


266 


1 


495 


2 


864 


3 


471 


2 


678 


3 


302 


TiOi 


mag 





022 





022 


0.016 


0.020 





024 





025 





022 





032 





031 





016 





023 





032 





027 





031 


Ti0 2 


mag 





053 





053 


0.041 


0.051 





056 





056 





052 





073 





070 





042 





056 





074 





062 





070 


TT £" 


X 

A 


2 


656 


2 


405 


3.501 


2.829 


3 


544 


2 


447 


2 


608 


-1 


280 


-0 


514 


3 


409 


2 


176 


-0 


668 


1 


180 


-0 


445 


H7A 


X 

A 


-0 


344 


-0 


822 


0.128 


-0.668 





334 


-1 


528 


-0 


705 


-4 


969 


-4 


202 





026 


-1 


476 


-4 


852 


-2 


984 


-4 


172 


HJf 


x 

A 


2 


676 


2 


265 


2.411 


2.261 


3 


096 


2 


084 


2 


287 





625 





934 


2 


373 


1 


973 





823 


1 


501 





944 


H7F 


x 

A 


1 


562 





643 


0.225 


0.347 


1 


653 





101 





559 


-1 


092 


-0 


865 





211 





035 


-1 


279 


-0 


531 


-0 


867 


Can A8498 


x 

A 





794 





854 


0.757 


0.820 





806 





842 





842 


1 


150 


1 


155 





775 


1 


001 


1 


122 





880 


1 


147 


Can A8542 


x 

A 


2 


723 


2 


741 


2.488 


2.673 


2 


615 


2 


710 


2 


720 


3 


334 


3 


342 


2 


514 


2 


939 


3 


259 


2 


757 


3 


323 


Can A8662 


2 

A 


2 


719 


2 


802 


2.683 


2.777 


2 


696 


2 


779 


2 


798 


3 


061 


3 


082 


2 


691 


2 


940 


2 


982 


2 


759 


3 


062 


Mgi A8807 


A 





448 





473 


0.343 


0.415 





382 





449 





460 





372 





382 





359 





343 





339 





465 





371 


CNB 


mag 





115 





080 


0.028 


0.060 





087 





051 





071 





202 





166 





030 





073 





161 





080 





163 


H + K 


mag 





333 





248 


0.183 


0.216 





258 





229 





235 





373 





336 





187 





239 





329 





263 





331 


Cai 


mag 





007 





009 


0.006 


0.007 





006 





005 





008 





015 





014 





007 





009 





013 





007 





014 


G 


mag 





173 





141 


0.073 


0.110 





134 





134 





128 





232 





211 





077 





132 





207 





168 





206 


Hb 


mag 





104 





046 


-0.036 


0.007 





073 





013 





032 





073 





064 


-0 


032 





020 





051 





026 





062 


MgG 


mag 





067 





064 


0.040 


0.060 





071 





067 





062 





102 





096 





042 





069 





107 





079 





096 


T\ TTT 

MH 


mag 





034 





036 


0.018 


0.024 





023 


-0 


028 





032 





060 





055 





020 





031 





051 


-0 


014 





055 


FC 


mag 





062 





056 


0.037 


0.049 





051 





052 





053 





073 





070 





038 





054 





069 





060 





070 


NaD 


mag 





057 





044 


0.014 


0.044 





053 





039 





043 





080 





075 





019 





057 





080 





053 





076 


D4000 




1 


785 


1 


572 


1.348 


1.481 


1 


617 


1 


543 


1 


533 


2 


070 


1 


938 


1 


364 


1 


576 


1 


971 


1 


677 


1 


924 


D4000„ 




1 


578 


1 


401 


1.223 


1.320 


1 


427 


1 


365 


1 


366 


1 


868 


1 


737 


1 


234 


1 


397 


1 


758 


1 


482 


1 


724 



Table 7: Values of the continuum indices for every galaxy class, continued from Tab. [6] 







l \ 1.1/ i 


111.!/ Z 


1 1 1. 1 / o 


RPD A 

1 1 1. 1 / ' 1 


1 1 1. 1 ' o 


□r 1 


or z 


or o 




or o 


CNj 


mag 


-0.002 


0.032 


0.063 


0.079 


0.066 


-0.047 


-0.069 


-0.064 


-0.054 


-0.047 


CN 2 


mag 


0.043 


0.078 


0.110 


0.124 


0.110 


-0.045 


-0.049 


-0.039 


-0.027 


-0.019 


Ca4227 


A 


1.038 


1.209 


1.304 


1.301 


1.286 


0.205 


0.341 


0.429 


0.539 


0.622 


G4300 


A 


4.289 


5.175 


5.563 


5.580 


5.465 


0.092 


0.597 


1.130 


1.659 


2.095 


Fe4383 


A 


4.177 


4.809 


5.087 


5.108 


5.046 


1.123 


1.194 


1.576 


2.064 


2.461 


Ca4455 


A 


1.418 


1.575 


1.604 


1.574 


1.536 


0.096 


0.361 


0.576 


0.735 


0.854 


Fe4531 


A 


3.190 


3.426 


3.499 


3.483 


3.430 


1.183 


1.625 


1.890 


2.122 


2.318 


C4668 


A 


5.160 


6.034 


6.615 


6.927 


6.669 


0.683 


1.198 


1.751 


2.337 


2.773 


H/3 


A 


1.407 


1.734 


1.597 


1.434 


1.409 


-12.371 


-6.644 


-4.584 


-3.640 


-3.117 


Fe5015 


A 


4.305 


4.728 


4.825 


4.741 


4.567 


-16.711 


-3.991 


0.305 


1.743 


2.286 


Mgi 


mag 


0.072 


0.089 


0.109 


0.120 


0.113 


0.019 


0.019 


0.024 


0.032 


0.039 


Mg 2 


mag 


0.186 


0.218 


0.251 


0.268 


0.256 


0.054 


0.068 


0.083 


0.101 


0.115 


Mg b 


A 


3.210 


3.626 


4.039 


4.232 


4.085 


1.301 


1.520 


1.763 


2.080 


2.333 


Fe5270 


x 

A 


2.544 


2.690 


2.706 


2.690 


2.669 


0.800 


1.229 


1.495 


1.710 


1.842 


Fe5335 


A 


2.370 


2.447 


2.380 


2.316 


2.320 


0.963 


1.295 


1.525 


1.708 


1.802 


Fe5406 


A 


1.525 


1.594 


1.582 


1.543 


1.525 


0.594 


0.800 


0.933 


1.054 


1.130 


Fe5709 


A 


0.889 


0.923 


0.901 


0.878 


0.897 


0.376 


0.505 


0.578 


0.654 


0.705 


Fe5782 


A 


0.770 


0.780 


0.774 


0.776 


0.785 


0.303 


0.412 


0.523 


0.631 


0.716 


Na D 


A 


2.732 


3.077 


3.575 


4.066 


4.481 


0.015 


0.633 


1.162 


1.898 


2.826 


TiOi 


mag 


0.026 


0.030 


0.033 


0.035 


0.034 


0.012 


0.014 


0.016 


0.017 


0.018 


Ti0 2 


mag 


0.062 


0.069 


0.076 


0.080 


0.078 


0.032 


0.034 


0.040 


0.046 


0.050 


tt r 


A 


0.762 


-0.698 


-1.724 


-1.993 


-1.728 


3.012 


3.859 


3.591 


3.057 


2.602 




A 


-2.713 


-4.343 


-5.508 


-5.856 


-5.574 


-1.284 


0.332 


0.260 


-0.311 


-0.868 


RSf 


A 


1.568 


0.893 


0.379 


0.222 


0.341 


1.653 


2.489 


2.463 


2.217 


1.986 




A 


-0.081 


-0.834 


-1.457 


-1.655 


-1.514 


-1.718 


-0.004 


0.301 


0.169 


-0.037 


Can A8498 


A 


1.078 


1.167 


1.144 


1.127 


1.187 


0.470 


0.561 


0.752 


0.807 


0.814 


Can A8542 


A 


3.231 


3.370 


3.318 


3.272 


3.445 


1.967 


2.166 


2.476 


2.581 


2.556 


Can A8662 


A 


3.074 


3.114 


3.010 


2.961 


3.182 


2.329 


2.476 


2.682 


2.728 


2.676 


Mgi A8807 


A 


0.410 


0.392 


0.365 


0.342 


0.478 


0.101 


0.200 


0.342 


0.408 


0.421 


CNB 


mag 


0.125 


0.175 


0.213 


0.229 


0.214 


-0.006 


0.021 


0.027 


0.037 


0.046 


H + K 


mag 


0.310 


0.352 


0.372 


0.375 


0.368 


0.113 


0.163 


0.182 


0.198 


0.218 


Cai 


mag 


0.012 


0.013 


0.016 


0.016 


0.016 


0.002 


0.005 


0.006 


0.007 


0.008 


G 


mag 


0.193 


0.224 


0.232 


0.227 


0.224 


0.030 


0.052 


0.071 


0.090 


0.106 


Hb 


mag 


0.067 


0.071 


0.067 


0.065 


0.065 


-0.166 


-0.070 


-0.036 


-0.020 


-0.010 


MgG 


mag 


0.080 


0.093 


0.108 


0.116 


0.110 


0.026 


0.033 


0.039 


0.047 


0.054 


MH 


mag 


0.044 


0.054 


0.063 


0.068 


0.065 


-0.064 


0.000 


0.020 


0.028 


0.033 


FC 


mag 


0.066 


0.071 


0.074 


0.075 


0.074 


0.017 


0.029 


0.037 


0.043 


0.046 


NaD 


mag 


0.060 


0.071 


0.085 


0.098 


0.105 


-0.031 


-0.006 


0.011 


0.029 


0.050 


D4000 




1.814 


1.993 


2.105 


2.146 


2.145 


1.161 


1.273 


1.342 


1.415 


1.514 


D4000„ 




1.611 


1.787 


1.903 


1.941 


1.919 


1.081 


1.166 


1.217 


1.271 


1.334 



Table 8: Value of the continuum absorption indices for the refined colour classes. 





iform [Gyr] 


7 [Gyr] 


*m [Gyr] 


t r [Gyr] 




TV 


logSSFR 


log 


M 
L 




median 


mode 


median 


mode 


median 


mode 


median 


mode 


median 


mode 


median 


mode 


median 


mode 


median 


mode 


p 


10.33 


11.42 


0.68 


0.89 


6.42 


6.42 


5.31 


4.57 


1.09 


0.89 


0.80 


0.06 


-12.11 


-11.40 


0.49 


0.49 


Ha 


10.33 


11.69 


0.66 


0.97 


6.42 


6.67 


5.55 


4.57 


1.21 


1.09 


1.16 


0.67 


-11.93 


-11.40 


0.57 


0.57 


SF 


8.69 


10.88 


0.50 


0.81 


5.17 


4.17 


4.08 


3.10 


1.09 


1.71 


1.90 


1.53 


-11.05 


-10.70 


0.57 


0.57 


A+Hn 


9.78 


12.78 


0.60 


0.85 


5.92 


5.42 


5.06 


4.57 


1.24 


1.71 


1.53 


1.28 


-11.58 


-11.05 


0.57 


0.57 


L 


10.33 


12.78 


0.66 


0.85 


6.67 


6.67 


5.55 


4.57 


1.28 


1.55 


1.28 


1.28 


-11.93 


-11.23 


0.57 


0.57 


S 


9.78 


12.78 


0.62 


0.77 


6.17 


5.42 


5.06 


4.57 


1.17 


1.71 


1.41 


1.28 


-11.58 


-11.05 


0.57 


0.57 


all 


10.33 


11.97 


0.68 


0.97 


6.42 


6.67 


5.55 


4.57 


1.21 


1.09 


1.04 


0.67 


-12.11 


-11.40 


0.57 


0.57 


P 


8.97 


10.88 


0.62 


0.83 


5.17 


4.67 


3.83 


3.34 


0.89 


0.62 


0.80 


0.06 


-11.40 


-11.05 


0.31 


0.31 


Ha 


9.24 


11.97 


0.58 


0.77 


5.42 


4.67 


4.32 


3.83 


1.01 


0.55 


1.28 


1.16 


-11.23 


-10.88 


0.40 


0.40 


SF 


7.33 


4.87 


0.44 


0.13 


4.17 


3.92 


3.10 


2.36 


1.01 


0.55 


1.65 


1.28 


-10.53 


-10.35 


0.40 


0.31 


A+Hn 


8.15 


6.24 


0.50 


0.30 


4.92 


4.67 


3.83 


3.34 


0.97 


0.55 


1.53 


1.16 


-10.88 


-10.53 


0.40 


0.40 


L 


9.51 


11.97 


0.58 


0.77 


5.67 


4.67 


4.57 


3.83 


1.01 


0.74 


1.28 


1.16 


-11.40 


-11.05 


0.49 


0.49 


S 


8.69 


6.24 


0.54 


0.77 


5.17 


4.67 


4.08 


3.83 


0.97 


0.55 


1.41 


1.16 


-11.05 


-10.70 


0.40 


0.40 


all 


9.24 


11.97 


0.58 


0.77 


5.42 


4.67 


4.32 


3.83 


1.01 


0.55 


1.16 


1.16 


-11.40 


-10.88 


0.40 


0.40 


P 


7.87 


6.78 


0.54 


0.91 


4.17 


3.42 


2.85 


2.12 


0.62 


0.20 


0.92 


0.06 


-11.05 


-10.53 


0.14 


0.14 


Ha 


7.33 


4.87 


0.44 


0.32 


4.17 


3.92 


3.10 


2.36 


0.86 


0.16 


1.28 


1.16 


-10.53 


-10.35 


0.31 


0.31 


SF 


4.33 


3.24 


0.34 


0.01 


2.42 


1.92 


1.63 


1.38 


0.93 


0.12 


1.28 


1.16 


-9.82 


-9.82 


-0.03 


-0.03 


A+Hn 


6.51 


4.87 


0.38 


0.03 


3.92 


2.92 


2.85 


2.12 


0.93 


0.16 


1.41 


1.28 


-10.35 


-10.35 


0.23 


0.23 


L 


7.33 


5.42 


0.46 


0.32 


4.17 


3.42 


3.10 


2.85 


0.82 


0.16 


1.16 


1.16 


-10.70 


-10.53 


0.23 


0.23 


S 


7.06 


4.87 


0.42 


0.13 


4.17 


3.42 


3.10 


2.36 


0.89 


0.16 


1.28 


1.16 


-10.53 


-10.35 


0.31 


0.31 


all 


6.78 


4.87 


0.40 


0.32 


4.17 


3.42 


3.10 


2.36 


0.89 


0.16 


1.28 


1.28 


-10.53 


-10.35 


0.23 


0.23 


P 


10.06 


11.15 


0.66 


0.89 


6.17 


6.42 


5.06 


4.57 


1.09 


0.89 


0.80 


0.06 


-11.93 


-11.23 


0.49 


0.49 


Ha 


9.78 


11.15 


0.62 


0.77 


5.92 


4.67 


4.82 


3.83 


1.13 


1.09 


1.16 


1.16 


-11.58 


-11.05 


0.49 


0.49 


SF 


4.87 


3.24 


0.34 


0.01 


2.67 


2.17 


1.87 


1.63 


0.93 


0.12 


1.28 


1.16 


-10.00 


-9.82 


0.06 


0.06 


A+Hn 


7.87 


5.69 


0.48 


0.30 


4.67 


3.92 


3.59 


2.85 


0.97 


0.55 


1.41 


1.16 


-10.70 


-10.53 


0.40 


0.40 


L 


9.78 


12.78 


0.62 


0.85 


6.17 


6.67 


5.06 


4.57 


1.24 


1.52 


1.41 


1.28 


-11.76 


-11.05 


0.57 


0.57 


S 


8.69 


6.24 


0.54 


0.30 


5.17 


4.67 


4.08 


3.83 


1.01 


0.55 


1.41 


1.16 


-11.05 


-10.70 


0.40 


0.40 


all 


9.78 


11.15 


0.62 


0.77 


5.92 


4.67 


4.82 


3.83 


1.13 


1.09 


1.16 


1.16 


-11.58 


-11.05 


0.49 


0.49 


RED 1 


9.24 


10.88 


0.58 


0.77 


5.42 


4.67 


4.32 


3.83 


0.97 


0.55 


1.16 


1.16 


-11.23 


-10.88 


0.40 


0.40 


RED 2 


10.06 


11.15 


0.64 


0.77 


5.92 


4.67 


4.82 


3.83 


1.09 


0.89 


1.04 


0.67 


-11.76 


-11.23 


0.49 


0.49 


RED 3 


10.33 


11.97 


0.68 


0.97 


6.42 


6.67 


5.55 


4.57 


1.17 


1.09 


1.04 


0.67 


-12.11 


-11.40 


0.49 


0.49 


RED 4 


10.60 


12.78 


0.70 


0.83 


6.92 


6.67 


6.04 


5.55 


1.28 


1.17 


1.16 


1.28 


-12.11 


-11.76 


0.57 


0.57 


RED 5 


10.60 


12.78 


0.68 


0.85 


7.17 


7.92 


6.29 


5.55 


1.48 


1.94 


1.53 


1.53 


-12.11 


-11.76 


0.74 


0.74 


SF 1 


1.33 


1.06 


0.44 


0.56 


0.67 


0.42 


0.40 


0.40 


0.86 


0.78 


1.16 


0.92 


-8.94 


-8.94 


-0.63 


-0.63 


SF 2 


2.69 


1.87 


0.36 


0.01 


1.42 


0.92 


0.89 


0.89 


0.93 


0.12 


1.16 


0.80 


-9.47 


-9.47 


-0.28 


-0.37 


SF 3 


4.06 


2.96 


0.34 


0.01 


2.42 


1.92 


1.63 


1.38 


0.97 


0.12 


1.28 


1.16 


-9.82 


-9.82 


-0.03 


-0.03 


SF 4 


5.69 


3.51 


0.36 


0.13 


3.42 


2.67 


2.36 


1.63 


1.01 


0.78 


1.41 


1.16 


-10.17 


-10.00 


0.23 


0.23 


SF 5 


7.60 


4.87 


0.46 


0.30 


4.67 


3.92 


3.34 


2.36 


1.01 


0.55 


1.65 


1.53 


-10.70 


-10.35 


0.40 


0.40 



Table 9: Modes and medians of the probability distributions of the physical parameters of the composites resulting from the Bayesian continuum 
fitting technique. ti mra : age of the continuous model component; 7: time-scale parameter of the star formation rate; tu'- average age of the stellar 
population, weighted by the stellar mass; t r : average age of the stellar population, weighted by the r-band magnitude; Z: metallicity; ry: optical 
depth of the dust component; SSFR: specific star formation rate; M/L: mass-to-light ratio. 



physical interpretation. We have calculated the most im- 
portant spectral feature indicators of the composites and 
shown that the parameters of the spectra span almost 
as wide a range as the spectral parameters of individ- 
ual galaxies do. This makes the atlas a representative 
collection of the numerous types of galaxies observed by 
SDSS. 

The composite spectrum atlas is available online at 
http : //www . vo . elte . hu/compositeatlas. 
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